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We propose a two-pulse excitation to efficiently photoinduced a novel sp 3 -bonded nano-domain 
with a frozen shear in a graphite crystal. This sp 3 structure is known to be well stabilized by shear 
displacement between neighboring two graphite layers. This shear motion is generated only as a 
transient and unfrozen one by the first visible laser pulse, shone over the graphite crystal. While, 
the second pulse is proved to freeze it before it disappears, resulting in an efficient interlayer a (sp 3 ) 
bond formation. 

PACS numbers: 64.70. Nd,61. 48. Gh,31.15.xv,05.40.-a 



Plenty of experimental and theoretical studies for pho- 
toinduced structural phase transitions (PSPT) have re- 
vealed the properties of multi-stability hidden in var- 
ious materials [l|-l5|. Carbon is one of such materials 
and exhibits various condensed phases with sp 2 and sp 3 
structures As well known, graphite, which is the 
most stable phase of condensed carbon systems, has a 
layer structure characterized by an sp 2 -bonded network. 

Recently, using scanning tunneling microscopy (STM), 
Kanasaki et al. have discovered that the irradiation of 
visible photons onto the graphite induces novel nanoscale 
sp 2 -)> sp 3 collective conversion^, [Tcj. The STM im- 
age shows the nanoscale domain in which 4 carbons in 
each six-membered ring rise up from the layer and resid- 
ual 2 sink down and close to the second layer, as shown 
in FigQJa). The sinking down carbons form interlayer 
a bonds between the first and the second layers like 
FigOJb). Due to the AB stacking of the graphite, there 
are two types of interlayer bonded carbon pairs, namely, 
a and f3 pairs. In order to stabilize these interlayer 
bonds, shear displacement occurs between the two lay- 
ers. We briefly recapitulate the essential aspects of this 
experiment. The exciting light with an energy of 1.57 eV 
should be polarized perpendicular to the graphite layers 
and should be a femtosecond (fs) pulse. A polarization 
parallel to the layers or a picosecond (ps) pulse gives 
almost no contribution. Thus, this nonequilibrium phe- 
nomena requires an interlayer charge-transfer excitation 
and a transient generation of an excited electronic wave 
packet. Moreover, the process is nonlinear but at a lower 
level than a ten-photon process. The resultant domain is 
stable for several days at room temperature. 

This phenomenon is closely related with graphite- 
diamond conversion, which has a long history of consid- 
erable experimental and theoretical researches [ii [Ti1 - [Tt| . 
The diamond is about 0.02 eV/carbon higher than the 
graphite [III, and the adiabatic energy barrier between 
them is about 0.3-0.4 eV/carbon[l, El, E3- Therefore, 
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FIG. 1: The structure of diaphite. (a) The top-down view of 
the first layer. The blue spheres indicate rising up carbons 
from the layer and the red ones are sinking down carbons, 
(b) The side view of the diaphite. When the first and the 
second layers just approach each other, the interlayer bonds 
in f3 pairs are unstable (upper illustration). Therefore, the 
shear displacement occurs so that the diaphite structure are 
stabilized (lower illustration). 



to achieve uniform transformation from the graphite to 
the diamond, a macroscopic order of energy is required. 
For this reason, the conventional graphite-diamond con- 
version is well known to occur only under high temper- 
ature and pressure (3000 °C JL _ 15 GP)[II EH, or by the 
irradiation of strong X rays 

Thus, the conventional conversion is quite different 
from the present photoinduced phenomenon. In the 
latter case, only a microscopic energy of several pho- 
tons, enough to nucleate a minimum kernel of the new 
phase, is given to a local and limited area of the crys- 



tal. This minimum kernel will proliferate stepwise, ac- 
cording to the hidden stability intrinsic to the graphite. 
The photoinduced nanoscale domain shown in Fig. [T]is 
not the conventional diamond, but an intermediate state 
between the graphite and the diamond, which is called 
"diaphite.fi, [l0|" 

From the above STM experiment, we can think of the 
following scenario for the present process [l8| . When a fs- 
laser pulse is irradiated onto graphite layers perpendicu- 
larly, an electron-hole pair spanning two layers is gener- 
ated. This electron-hole pair mainly dissipates into the 
semimetallic continuum of the graphite as plus and minus 
carriers due to the good conductivity of graphite. How- 
ever, by a small but finite probability, this electron- hole 
pair is expected to be bounded with each other through 
the interlayer Coulomb attraction. This exciton-like state 
self-localizes at a certain point of the layers by contract- 
ing the interlayer distance only around it. As the lo- 
cal contraction of the interlayer distance proceeds, an 
interlayer a bond is formed. Through further pulse ex- 
citations, an increasing amount of interlayer a bonds are 
formed stepwise, and then the diaphite structure is ex- 
pected to appear macroscopically. 

We have already clarified the adiabatic property of 
the phase transition by calculating the adiabatic path 
from the graphite to the small diaphite domain, by 
means of the LDA[l9[ and by the semiempirical Bren- 
ner's theory [20I [2l[ which can describe various carbon 
cluster systems. Next, we have studied the early stage 
dynamics of the interlayer a bond formation along with 
the above scenario [HQ. As a result, we found that the 
self-localization of an electron-hole pair spanning two lay- 
ers occurs within a few femtoseconds at the probability 
of about 2% when the electron- hole pair is excited as a 
transient state with the energy of 3.3 ± 1.8 eV. After the 
self-localization, by using classical molecular dynamics 
(MD) with the Brenner's potential, we obtained that the 
subsequent bond formation is achieved within about 0.5 
ps when the excitation energy is more than 4.5 eV, corre- 
sponding to about three visible photons with the energy 
of 1.57 eV. 

Moreover, we investigated the effects of random multi- 
site excitation by using the MD[23[. In this calculation, 
we select the excited sites randomly according to exci- 
tation density, and give random energy to each excited 
site according to the Gauss distribution with the mean 
energy of 3.3 eV and the width of 1.8 eV. Consequently, 
due to the cooperative phenomena between excited sites, 
the number of formed interlayer bonds increases nonlin- 
early with respect to the number of excited sites. The 
nonlinearity shows 1.7 power of the number of initially 
excited sites, indicating three or four photons are effec- 
tively involved in this interlayer bond formation. This 
theoretical result is also confirmed qualitatively by the 
experiment [q[ . 

However, the structure of the obtained domain with in- 
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FIG. 2: (a) The snapshots of the MD calculation in the case 
of 6% excitation at t=0, 0.4, 0.8 and 1.4 ps. These are the 
top-down view of graphite layers. The system consists of two 
layers with AB stacking and the initial interlayer distance 3.35 
A. Each layer includes 6240 carbons in about 128 A x 128 A, 
wherein a periodic boundary condition is imposed. The shear 
between the upper layer and the lower one and its direction 
are represented by the circular color palette at the top. (b) 
The average shear displacement as a function of time. It is 
averaged over all carbon sites and over 50 calculations. 



terlayer a bonds is quite different from the diaphite. In 
our calculation, all the interlayer bonds have been formed 
only in a pairs (we call them a bonds), not in f3 pairs 
(P bond) at all. As shown in Fig. Hfb), the diaphite 
includes a and j3 bonds alternatively. The j3 bonds are 
quite unstable due to the AB stacking, so that they have 
not been formed in our calculation. As mentioned above, 
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the ft bonds are stabilized by the shear displacement be- 
tween layers like the lower structure of Fig. QJb). The 
existence of the shear is particularly important for the 
formation of the diaphite domain, as already discussed 
in our previous works It has been experimen- 

tally confirmed before [9, 24] that the shear displacement 
is generated under fs-laser irradiation. However, it has 
not yet been clearly known how the shear occurs and /? 
bonds are successfully formed. 

We discussed in our previous paper that when the 
graphite is excited only at t = the shear once appears 
due to the excitation of /3 pairs but disappears with time, 
indicating that no /3 bond is formed [23j. Fig. [2fa) shows 
the snapshots of the MD calculation in the random mul- 
tisite excitation, representing the shear displacement of 
the upper layer relative to the lower one together with 
its direction by a circular color palette at the top. The 
excitation density is 6 %, corresponding to 374 carbon 
pairs' excitation. In this calculation, assuming only the 
x-polarized light irradiation in xy plane, all the excited 
/? pairs have the shift only in the x direction. 

As seen from these snapshots, the shear in +x direc- 
tion occurs in all the area once but afterward it dis- 
appears. Fig. [2fb) is the average shear displacement 

A£ = iV-jEtr [^eEf=r A*,], as a function of 
time. Here, N ca i is the number of calculations (=50), 
N s ue the number of carbon sites in each layer, and 
Axj = x u j — x\j, where # u (i),j is the x coordinate of 
the j th site in the upper (lower) layer. In the 6% excita- 
tion case, the average shear grows maximum at about 0.8 
ps. Thus, the shear is proved to be only a transient and 
unfrozen one, unless the /? bonds are formed. However, 
since the shear makes f3 bonds stable, they are expected 
to be formed more efficiently if the graphite were excited 
again before this transient shear disappears. In order 
to confirm it, we here investigate the effect of two-pulse 
excitation, using the MD. 

The calculation procedure and the system are same as 
in the random multisite excitation case. After a certain 
period of time At passes from the first excitation, we 
give additive excitation to the system. The excitation 
density of the second excitation is same as that of the 
first excitation, where we consider only 6% excitation 
which is enough to induce cooperative phenomena in the 
interlayer bond formation [23|. We count out the number 
of a and f3 bonds at 1.0 ps after the second excitation. 
By performing the calculations for various values of the 
time interval A£, we investigate the At dependence on 
the number of these two types of bonds. 

As a example, we show the snapshots of the MD calcu- 
lation for the case of At = 0.4 ps in Fig. [3l These snap- 
shots represent the interlayer distance by a color palette 
at the top, where the a and f3 bonds are colored by red 
and blue, respectively. Only the snapshot at the bottom 
right represents the shear displacement at 1.4 ps. As seen 
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FIG. 3: The snapshots of the MD calculation at t=0.4, 0.5, 
0.8, 1.1 and 1.4 ps in the case where the additive excitation 
with 6% density is given to the system at £=0.4 ps. The 
distance between the two layers is represented by the color 
palette at the top. The a and f3 bonds are colored by red and 
blue, respectively. Only the snapshot at the bottom right 
shows the shear at £=1.4 ps. 



from this figure, although there is almost no j3 bond until 
0.4 ps, we can see a lot of j3 bonds after the second exci- 
tation at 0.4 ps. These f3 bonds tend to be formed near a 
bonds, since the interlayer distance is contracting there. 
Even at 1.4 ps, some of them remain together with a 
bonds, resulting in the small domains in which the inter- 
layer distance is contracted. Moreover, as shown in the 
snapshot of the shear at the bottom right, the color of 
the domain including f3 bonds remains red. This means 
that the shear does not disappear with time, and remains 
around the (3 bonds. In other words, the transient shear 
is frozen by the second excitation. Thus, we confirmed 
that /3 bonds can be formed due to the two-pulse excita- 
tion. 
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layers by means of the MD calculation. By the fs-laser 
excitation of graphite, the interlayer contraction and in- 
plane shear displacement are generated only as coherent 
and transient lattice motions. These lattice motions are 
getting frozen as interlayer sp 3 bonds increase stepwise 
by the additive pulse excitations. We have thus proposed 
the two-pulse excitation to achieve the efficient transition 
from the graphite to the diaphite. 
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FIG. 4: The number of interlayer bonds counted at 1.0 ps 
after the second excitation, as a function of the time interval 
At between the first excitation and the second one. The red 
and blue solid lines indicate the number of a and f3 bonds, 
respectively, and the green dashed one is the sum of them. 
The values at At=0 mean no additive excitation, correspond- 
ing to the 12% excitation only at once. The black rhombus 
at At=0 corresponds to the 6% excitation only at once. 



Fig. H] shows the At dependence on the number of 
interlayer bonds, where the red and blue solid lines in- 
dicate the number of a and f3 bonds, respectively, and 
the green dashed one is the sum of them. The values 
at At=0 mean no additive excitation, corresponding to 
the 12% excitation only at once. For reference, the black 
rhombus at At=0 indicates the number of a bonds for 
the 6% excitation at once. As seen from this figure, the 
number of a bonds monotonically decreases and settles 
down to a constant at more than At =0.9 ps, while /3 
bonds increase until At =0.5 ps and afterward decreases 
gradually. In this calculation, the f3 bonds are formed 
most efficiently around At =0.3 - 0.5 ps. In this time re- 
gion, the average shear is most rapidly increasing as seen 
from Fig. [2](b). Although the number of /? bonds are 
still quite less than that of a bonds even in this region, 
it is essential that j3 bonds can be formed by the addi- 
tive excitation. Because, if the j3 bonds are once formed, 
the shear also remains stable around them. Therefore, 
through further excitation, new /? bonds will be formed 
more easily around the existing shear. Thus, f3 bonds will 
proliferate comparable to a bonds will proliferate. As the 
number of these two types of bonds increases enough, the 
diaphite will be achieved. 

In conclusion, we have theoretically studied the 
nonequilibrium formation process of sp 3 -bonded carbon 
nanodomains with shear displacement between graphite 
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